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Abstract. We analyze the Jordan-Brans-Dicke model (JBD) of gravity, where 
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deviations from General Relativity (GR) are described by a scalar field non-minimally 



O ■ coupled to the graviton. The theory is characterized by a constant coupling parameter, 

w , wjbd; GR is recovered in the limit cjjbd oo. In such theories, gravity modifications 

manifest at early times, so that one cannot rely on the usual approach of looking 
for inconsistencies in the expansion history and perturbations growth in order to 



discriminate between JBD and GR. However, we show that a similar technique can 
■ be successfully applied to early and late times observables instead. Cosmological 

parameters inferred extrapolating early-time observations to the present will match 
those recovered from direct late-time observations only if the correct gravity theory 
^) ■ is used. We use the primary CMB, as will be seen by the Planck satellite, as the 

early-time observable; and forthcoming and planned Supernovae, Baryonic Acoustic 
Oscillations and Weak Lensing experiments as late-time observables. We find that 
detection of values of wjbd as large as 500 and 1000 is within reach of the upcoming 
(2010) and next-generation (2020) experiments, respectively. 



1. Introduction 

Theoretical modifications to General Relativity (GR) were in many cases formulated 
much before the discovery of cosmic acceleration, often within the more general quest 
of the search for a parent theory of GR which could lead to the Grand Unification of 
the four fundamental forces (e.g. (jl|; 0) and references therein). In the last few years, 
however, possible applications of modified gravity theories as an explanation for cosmic 
acceleration have generated a renewed interest in the cosmological community towards 
such models, and many new ones have been formulated (e.g. fi B B B Q B B B Q 
and references therein). The most compelling shortcoming of the standard Cold Dark 
Matter model is the observed mismatch in the right and left hand side of Einstein 
equations, which can be mostly simply accounted for through a Cosmological Constant 
term. However, the form of such equations is derived from the Lagrangian of GR, so 
that an alternative explanation is that the latter breaks down on cosmological scales. 
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In this paper we analyze possible observational sig natures of a class of modified GR 



models known as scalar-tensor theories (e.g. (Il2l ; Il3l ; Il4l)). where the action of gravity is 
determined by a scalar field in addition to the metric tensor. We consider the simp lest 
example of scalar-tensor theory, the Jordan-Brans-Dicke model (JBD), (1 15c Il6l ); in 
this case, the coupling between gravity and the scalar field is described by a constant 
parameter, wjbd- The main motivation for this work is to provide a general approach to 
test GR on cosmological scales, as opposed to Solar System scale measurements; this is 
in principle a well distinct task from that of addressing the cosmic acceleration problem. 
In particular, as we will show later, in the JBD model the modifications of gravity are 
limited to early times, so that acceleration cannot be obtained as a result of the gravity 
modification alone. However, the theory predicts the introduction of a scalar field as 
a gravitational degree of freedom, and such scalar field behaves as minimally coupled 
at late times, so that it can play the role of Quintessence, given a suitable potential. 
In this sense, the JBD model can provide an explanation of the nature of Dark Energy 
(1171 ); it does not address the fine tuning problem, which is, however, common to most 
dark energy models. 

The strongest constraint to date on this model has been put on the Solar System 
scale: the present 2a limit from the Cassini spacecraft is cujbd > 40000 (1181 ). However, 
such constraint does not necessarily apply on distances much larger than those of the 
measurements, and epochs much different from the present. Local universe experiments 
only probe scales in gravitational equilibrium, where the background expansion of 
the Universe is negligible; they would not reveal spatial or time variation of the 
gravitational constant on larger scales (jl9t |20| ). Sensible limits on the value of wjbd 
which is representative of the whole Universe have to be inferred from observations 
on cosmologically relevant scales (12 ll ). The current limit on such quantity, obtained 
combining the extended WMAP 1st year data and the 2dF large-scale structure data, 
is cujbd > 120 at 95% confidence level (j22j). 

In the present paper we will show how constraints on the "cosmological" JBD 
parameter can be substantially improved using suitable combinations of next-generation 
experiments. We consider the CMB power spectra, coupled at a time with Supernova? 
type la (SNe), Baryon Acoustic Oscillations (BAO), and Weak Lensing (WL). The 
main difference with respect to earlier works is that we will not assume previous 
knowledge of the model. We propose a general method for discriminating between 
GR and JBD, which relies on the fact that the two models are different at early times 
and similar at late times. Therefore, observations at early times, if extrapolated to 
the present epoch, will be consistent with late times observations only if the correct 
theory of gravity is used. This is somehow close in spirit to the consistency checks 
between expansion history and perturbations growth often used in order to discriminate 
between dark energy and modifications to GR, in the case that they differ at late epochs 
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32j), and references therein). 



(e.g. m 

The method presented in this paper can be easily applied to any pair of models 
which agree at some epoch and disagree at some other epoch. This include not only 
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modifications to GR as opposed to GR (1331 ; 134 ). but also, for example, Ordinary 



(35 
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371 ) versus Early (J38 
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), or Extended (1421 ; |43| ; |44|), Quintessence 



models. 

The outline of the paper is as follows. In Sec. [2] we present the JBD model and its 
general phenomenology. In Sec. [3] we compare a dark-fluid GR model with the same 
expansion history and the JBD model, showing how the JBD field behaves in terms 
of equation of state and how it would change the interpretation of the cosmological 
parameters if Cosmology was assumed to be GR. In Sec. 0] we present our method in 
detail, describe the observables that we are going to use, and present results for two 
values of u jbd of interest. Finally, in Sec. Ewe summarize our results and discuss their 
implications. 



2. Phenomenology of the JBD model 

The JBD cosmological model was formulated in 1961 as as the first scalar-tensor theory 
of gravity (flol ). It only features one more degree of freedom with respect to GR, the 
JBD parameter cujbd, which is constant both in space and time. The Lagrangian of the 
JBD model reads 

£jbd = ($R - - V(&) + £ fluid (1) 

where R is the Ricci scalar, and £fl U id is the Lagrangian of the ordinary matter and 
radiation components, . Its equation of motion are 

H +iJ ¥--6~l v ¥j + ~¥' (2) 
* + aF *=a5^3^-^> (3) 

here H is the Hubble parameter a /a, where a is the scale factor, and dots denote 
derivatives with respect to proper time. 

GR is a particular case of the JBD theory, corresponding to cujbd = oo. For such value 
the solution of the above equations is <fr = cost = 1, and <&/<!> — ► steeply, so that 
additional terms in the Friedmann equation disappears and its ordinary gravity form is 



recovered (1451 ). 



2.1. Background dynamics 

To study the phenomenology of the JBD model and its differences with GR, we start 
by comparing cosmological observables for the two models assuming that they have the 
same cosmological parameters other than cujbd- The GR ACDM model is specified by 
6 parameters: o;& = Qf,h 2 , ujcdm = ^CDAih 2 , n s , r, A s , h; physical density of baryons, 
physical density of cold dark matter, primordial power spectrum spectra slope, optical 
depth to the last scattering surface, amplitude of the primordial perturbations and 
Hubble constant in units of 100 km/s/Mpc. When comparing to a JBD model, we 
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assign the GR ACDM model a formal value of ujjbd — oo. 
It is useful to re-define the JBD field in the following way: 

With such notation the field 4> is now a canonical scalar field, non-minimally coupled to 
the Ricci scalar, and the Lagrangian of the JBD model reads 

£jbd = \^ 2 R - \d,<pdy - V{(j>) + £ fluid . (5) 

For comparison, the Lagrangian for GR with a Cosmological Constant term is: 

£gr = Y^q R - A + ^fluid- (6) 

To isolate the effect of the gravity modifications alone, we set the potential to be 
constant, mimicking a A term, throughout all of our analysis. This corresponds to the 
"worst-case scenario" in terms of detectability, since in this case Eqs. (JH]) and (jSJ) only 
differ in the gravity sector, and late-time acceleration is driven by the same mechanism. 
The equations of motion for the field are given by 

II- O U " - / "M i 4 WjBD C7\ 

H + 2H - = -^j BD I ~ ) + 3^2~^ ( 7 ) 

2o;jbd P - 3p 

2CUJBD + 3 2 U 

Initial conditions for the evolution of the field can be set requiring that its present value 
reproduces the strength of gravity observed in Cavendish-type experiments in the local 




universe (1451): 

(9) 



,2 _ ^JBD 2 C^JBD + 4 

~ 2ttG 2cujbd + 3' 



The above equation is ensuring that at the present time the gravitational coupling 
£0 2 — > l/(8nG), so that GR is the late-time limit of theory for any value of c<Jjbd Ly- 
The other initial condition is 0beg = 0, since the trajectory of the field during the 



radiation era is constant, as known from analytical solutions in this regime (1461 ; 1471 ). 
With the notation £0 2 = F(<j>), we can now explicitly compare the evolution of the 
Hubble factor for the JBD model and for GR: 

' Luid + ^ (U 2 + a 2 V- 3HF) 



H 



2 

JBD 



3F 



8nG J 



q 8ttG A . , 

Hgr = ^-Pfluid + 3 ■ (10) 

The first deviation from GR is expressed by the 1/3F term, which acts like an effective 
time- varying gravitational constant G c s = 1/(8ttF). Its evolution is shown in the first 

J There is indeed a higher order correction to the quoted limit, known as Cavendish correction, 
expressed by the second term on the right hand side of Eq. ([9]). We take it into account in our 
numerical codes, but we will generally say, slightly inappropriately, that the theory recovers GR at the 
present time. 
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Figure 1. Left and right panel: evolution of the gravitational coupling 1/F, 
normalized to the present value, and of the squared Hubble factor, normalized to 
a GR model with the same cosmological parameters. We show curves for wjbd = 100 
(solid line), 500 (dashed line), 1000 (dotted line). Center panel: relative contributions 
of the four terms in the r.h.s. of Eq. (fTOjl : fluid part (solid line), kinetic term (dotted 
line), constant potential (dashed line), F term, plotted with opposite sign (dashed- 
dotted-line) ; trajectories are shown for a fiducial model with wjbd = 100. 

panel of Fig. [TJ for fiducial models with wjbd = 100, 500, 1000. We have normalized the 
curves to their present value, which also correspond to the standard GR gravitational 
constant since the theory approaches GR at late times. It is a monotonic, decreasing 
function of time, implying that the gravitational force was stronger in the past. Its 
trajectory stays close to the initial value in the radiation and early matter domination 
era: for such redshifts the difference between GR and JBD is largest. For larger values 
of 

a; jbd the curves are closer to the GR value and their time variation is less significant. 
The second difference is caused by the additional terms in the right hand side of Eq. 
( ITU]) . For the case ujbd = 100, we plot their relative contributions to the total energy 
density in the middle panel of Fig. [TJ The energy density is governed by the fluid-like 
components (matter and radiation) up to late times, when the effective cosmological 
constant (constant potential) takes over: this is similar to what happens in GR. The 
only difference is a slight shift due to the term 3HF/a 2 , which has a relative weight of 
a few per thousand up to one per cent for this value of u;j B d, while the kinetic term is 
at least one order of magnitude smaller. 

We conclude that the largest impact of these modifications of gravity is caused by the 
change in the effective gravitational constant, which is proportional to l/F(<f)), and most 
relevant at early times. 

Finally, the right panel of Fig. [ljwe show the ratio H] BD / Hq R ; the different evolution 
of H 2 (z) encloses both the two effects described above. 

3. Interpretation of the JBD field: early-time effects 

We now proceed to understand how evidence of non-GR gravity can be discovered, with- 
out assuming any a priori knowledge of the gravity theory. 

A key element in this respect, that has been extensively recognized in literature, is the 
fact that modifying gravity affects both the background evolution and the linear per- 
turbations growth, n particular, in GR the redshift dependence of the matter density 
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perturbation (5 m = Sp m /p m , where p m denotes the matter density) can be exactly pre- 
dicted for a given expansion history Thus in general, if one could measure H{z) and 
S m (z) at the same time and infinitely well, any modification of the underlying gravity 
theory would manifest as an incongruence between the prediction of GR and the actual 



measurement (|23|; |2J; [25|; \M \M l2Sl; |2S l30; |3l|; |4g|). 



However, modifications to the redshift evolution of 5 m in JBD result in maximum part 
from the different strength of the gravitational field, which is encoded, as seen in the 
previous Section, in the function 1/F. Such function is significantly different from the 
GR case only at high redshift, since by construction the model tends to GR at late 
times. One would need to accurately measure the matter density perturbations at early 
times in order to distinguish between the two cases. Unfortunately, at such epoch direct 
measurements of S m (z) are expected to be extremely challenging. (More details on the 
evolution of 5 m in the JBD model can be found in Appendix A). 

We conclude that fixing the expansion history for such models also corresponds to have 
very similar perturbations growth rate at redshifts relevant for structure formation, so 
that the idea of detecting deviations from GR by means of these two observables cannot 
be applied in the present case. 



4. Early times vs late times observables 

Although JBD modifications to GR do not manifest as incongruencies between the 
expansion history and the growth of perturbations, a similar approach can still be used. 
We propose to use early and late time observables. 

The key element of this method resides in the fact that by construction our JBD 
models resemble GR in the local universe. Thus, a standard GR, ACDM model will 
be a good fit to low redshift data and will recover the correct underlying cosmological 
parameters whether or not the "true" theory is JBD or GR. Conversely, observations at 
high redshift would be sensitive to the modified gravitational coupling, and give different 
best fits for the same parameters, if gravity is JBD {ujjbd < oo), but agree with the low 
redshift parameter fit if gravity is GR. 

Although we present and develop the method in the context of the JBD theory, the 
power of such method is that no previous knowledge of the model is assumed; the only 
requirement is to use different datasets separately according to their "early" or "late" 
time nature. Moreover, even if we use the JBD model as a working example, and will 
give limits on the JBD parameter, such inconsistencies are indeed a general indication 
of deviation from GR at some level. 

4-1. Getting quantitative 

To quantify how small deviations from GR can be observed through the method 
described above, we will start by assuming that the "true" Universe is described by 
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a JBD theory with a given value of cujbd- We will then generate a typical early-time 
observable, associate to it expected observational uncertainties, and fit it with a GR 
ACDM model. This gives an estimate of the best fit values of the recovered cosmological 
parameters, as well as of errorbars. For all the models which lie within la from the 
best fit model, we then make predictions for the value of late-time observables, again 
assuming GR. We then compute the late-time observable in the "true" JBD model and 
its associated expected observational uncertainties. We obtain two different confidence 
regions for the same observables: the first is predicted from an early-time observation 
and the (wrong) assumption of standard gravity, while the second is the confidence 
region allowed by observations. If the difference has enough statistical significance, 
it means that observations are able to discriminate between the two cases, JBD and 
Einsteinian gravity, for that particular value of u;jbd • 



4-1.1. Early time observable For the early-time observable we choose the Cosmic 
Microwave Background (CMB), as it can be measured by the Planck experiment (|49l ). 
We use the Planck specifications as described in the Planck Blue Book (l50l ). 
Let us note that, although not an object of study of the present work, another 
possible early time observable would be the the amount of baryons from Big Bang 
Nucleosynthesis, which is modified in JBD scenarios with respect to GR (e. g. (15 ll : I52l0 . 
The most stringent limit obtained on u;jbd from BBN constraints is o;jbd > 32 (1531 ) . 
confirmed by the analysis in (54). Such limit improves significantly using the more 
recent, tighter constraints on the abundance of primordial deuterium to hydrogen 
ratio from (55|; (HHI) . Assuming that the current upper limit on D/H can be cast as 
D/H < 4 x 10~ 5 , the limit on u;jbd tightens to cujbd > 90. 

We consider several representative values for u JBD : 100, 200, 500, 1000. This 
range is motivated as follows. As for the lower bound, (1221 ) showed that the fitted 
values of cosmological parameters for a JBD case with o;jbd = 70 are already ruled out 
by data. For the upper bound, a Fisher matrix analysis for a JBD model and the Planck 
experiment forecasts a 1 — a detection threshold of the u;jbd parameters between 1000 
and 2000 (1571 ) . While the Fisher matrix approach already assumes knowledge of the 
underlying model (and thus give more stringent constraints than model independent 
methods), the method presented here relies in the additional statistical power of late 
time observables. 

For each of the reference JBD models, we compute the "true" CMB temperature, 
EE polarization, and their cross-correlation power spectra with the DEfast code 



(58 



42; 



431 ) . originally based on the CMBfast package (1591 ) . We find the best fit GR 



ACDM parameters and their confidence regions by running Markov Chain Monte Carlo 
chains, using the COSMOMC code (60), which uses CAMB (joll ). One of our concerns 
was that small numerical differences between the two codes could give rise to systematics 
effects in the recovered parameters. In fact, even running the two codes with settings 
as close as possible, some numerical differences in the CMB power spectra remain. We 
find that such difference is independent of the cosmological model, and we correct for 
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it, as described in detail in Appendix B. In Tabled] we show the best fit cosmological 
parameters obtained fitting CMB data for various wjbd; we report uncertainties at the 
2cr level. For reference, we also show how well we recover the cosmological parameters 



Table 1. Best fit and 2a confidence levels of cosmological parameters from Planck 
CMB forecasts, for JBD models with uj jbd = 100, 200, 500, 1000, oo, fitted with a 
ACDM model. A = ln(10 10 A s ). 





input 
JBD model 


wjbd = 100 
ACDM fit 


wjbd = 200 
ACDM fit 


wjbd = 500 
ACDM fit 


WJBD = 1000 

ACDM fit 


WJBD = OO 

ACDM fit 


UJb 

^CDM 
n s 

T 

A 

x 2 


0.022 
0.1232 
0.95 
0.09 
3.1355 
0.72 
0.28 



n noi c 0.0218 
U.UZIO 0.0212 
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U.IZOUq.1219 

0-9295 lifi 
0.0873o ° 7 9 | 7 8 

o inn 3.1333 
o.riiz 3 0991 

66^0 °- 6798 
U.ODOU 0.6508 

0.3295 '.3084 

24.64 


0.02198:811 

o.i24i8: 1 1 2 2 7 oI 
0.9448 8:111 

o.o9oo8;81!l 

3.1349 iiflt 

0.6965 8:S 27 
o.3oo98:|| 

9.96 


n noo 0.0223 

u - uzz 0.0217 

0.12368i 2 2 oi 
0.94788:E1E| 
0-0905 « 
3.1376|:i||g 

0.7088 '.6940 

0.2899 8:1??? 

9.22 


n noo 0.0223 

U.UZZ 0.0217 
n iooc 0.1269 
U.1ZJO 0.1202 

n QM 1 0.9588 
U.yoil 0.9429 

0-0915 g-.oli 
3.i4ii |:lf 9 7 

n 71 Q1 0.7346 
U-'iyi 0.7054 
n 901 c; 0.2989 
U.ZOIO 0.2638 

9.06 







for a ACDM model, formally identified by cjjbd — oo. 

Some of the parameters are quite insensitive to the presence of the JBD field, and their 
predicted values coincide with those of the "true" JBD model within errorbars, even 
for values of c<Jjbd as small as 100. This is the case for r and the primordial amplitude 
of the perturbations, A s . Interestingly, the combination u m = f2 CDM /i 2 + Vt h h 2 is also 
recovered quite well. In fact, such combination is mainly constrained by the distance to 
last scattering, which is in large part integrated over look-back times where the field is 
unimportant. However, h and Q m are respectively under- and over-estimated. In fact, 
the size of the sound horizon at recombination is significantly smaller in a JBD scenario 
than in a ACDM (e.g. (1621)). since it is only affected by pre-recombination physics, 
where the relative weight of the field component in the total density is significant. As 
we have seen in Sec. 3, gravity in GR is weaker at early times than in JBD and the JBD 
phenomenology is in part mimicked by an additional matter field. This moves the peaks 
towards smaller scales; since the distance to last scattering is constrained, such effect 
can be only accounted for enhancing the matter component, and lowering the Hubble 
factor. 

Note that even for a value of o;jbd as large as 500, the shift in Q m from its true value 
is as large as 7%, and that in h is between 3 and 4 %; furthermore, their "true" values 
lie outside the 2 — a confidence contours obtained through Planck-quality CMB data. 
This is interesting because these parameters can be constrained well through late-time 
observables. 

Finally, there is an overall shift of power towards smaller scales, corresponding to an 
underestimation of the primordial power spectrum index, n s , also resulting from the 
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Figure 2. Left panel: Temperature power spectrum of the reference JBD model 
with a; jbd = 100 (solid line), compared to a GR ACDM model with same cosmological 
parameters (dashed line) , and to the reconstructed GR A CDM model with parameters 
of Table Q] (dotted line). Although the reconstructed model resembles the "true" JBD 
model much more closely, some differences are still present. Right panel: Ratio of 
the reconstructed temperature power spectrum and the input "true" JBD model, for 
wjbd = 100 (solid line), and uj jbd = 500 (dashed line). The residual difference seen 
for ojjbd = 100 is reflected by the the high value of \ 2 of the fit. 



change in the location of the peaks. 

In the last entry of the table we report the x 2 f° r the best fit. The value of x 2 =9.06 for 
the GR case cjjbd — oo arises from numerical differences between DEfast and CAMB. 
For reference, the value of the x 2 f° r the GR model if the same code, either DEfast or 
CAMB, is used both for generating the CMB spectra and for the fitting procedure is of 
order unity or smaller (see Appendix B for more details). 

In Fig. [2] we show the temperature power spectrum for the reference JBD model 
with cjjbd — 100, for a ACDM one with the same parameters, and for the best-fit 
ACDM, with the parameters of Table [U There is a residual disagreement between the 
spectrum obtained from the fitting procedure and the "true" JBD spectrum, reflected 
by the relatively poor value of the x 2 °f the fit. In particular, differences are seen in 
the ISW amplitude and in the height of the first two peaks. The first is generated by 
the clustering properties of the JBD field, which cannot be reproduced with a smooth 
component; the second comes from the slight underestimation of the baryon content. 
The different amount of ISW could in principle be used in order to detect modifications 
to GR (42), but its signal-to-noise level is expected to be below the cosmic variance for 
the target values of ^jbd considered in this work. 

Such differences becomes negligible as long as cujbd is larger than 500, as can be seen 
in the right panel of Fig. [2j We show the ratio of the temperature power spectra of 
the GR ACDM model coming from the fit and of the true JBD model input, for wjbd 
= 100, 500 respectively. We conclude that for values of co>jbd smaller than or equal 
to 500, high values of the x 2 are already a sign that the fit is done using the wrong 
model. More sophisticated techniques of Bayesian analysis may be enough to show that 



more cosmological parameters are needed, e.g. (1631 ; 164 ; l65l ; l66l ) . However, in the regime 



where wjbd is larger than 500, the fit is almost as good as the one for the GR case. For 
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Figure 3. Linear fit of 1/wjbd versus f2 m , h and x 2 as recovered from a fit to a GR 
ACDM model. Data points are shown for wjbd = 100, 200, 500, 1000 and oo; the last 
case corresponds to GR. 

such values, one can not rely on the CMB alone, coupled to goodness-of-fit techniques; 
the consistency between early and late time observables needs to be used. This paper 
focuses on this regime. 

One interesting feature of the observed shift in the cosmological parameters, as 
reconstructed by the "wrong" fit, is that it scales almost linearly with 1/cjjbd, at least 
for the most sensitive Q m and h. The value of the x 2 a ls° follows a similar pattern, but 
the fit is different for values of cojbd — and > 500, since the \ 2 is almost constant above 
this value. We show such linear fits, of the form of par fit = a + b x (1/u;jbd), m Fig- 
[3] for reference, we also plot the 2a limits coming from the fit. Coefficients of the fit 
are a=(0.294, 0.719, 0.385) and 6=(10.52,-10.75, 4286) for Q m , h and xVjbd < 500) 
respectively; for larger values of cujbd the \ 2 is fitted by a p i ateau = 8.96; & p i a teau = 450.3. 
This property is useful because it allows to quantify in a simple way the effect of the 
JBD field on the cosmological parameters. In particular, it can be used to estimate 
which target cujbd experiments could detect, starting from the attainable precision on 
the cosmological parameters. 

4-1-2. Late time observables We use the Monte Carlo simulations described above and 
select all models within 68.3% confidence level from the best fit ACDM model; for these 
models we extrapolate the values of late-time observables. The late time observables we 
consider are: the distance modulus as seen through Supernovae observations, y,(z)\ the 
angular diameter distance oIa{z) and the Hubble factor H(z), by means of the Baryon 
Acoustic Oscillations; and the power spectrum of convergence P K , with weak lensing 
surveys. We compute the same late-time observables for the JBD fiducial models and 
their respective observational errors for different experimental set ups. We concentrate 
here on the possibility of constraining values of cujbd > 500, for which, as we said, we 
cannot use information about the quality of the fit, postponing the analysis of lower 
values of co> JB d to future work (joTl ). we therefore consider, at first, datasets that will 
be available around 2010 ("2010"); and second, datasets likely to be available a decade 
later ("2020"). 

We evaluate observational errors, for each observable, as follows; specifics of the exper- 
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imental setups considered can be found in Tables [2J [31 lU 



SNe. For Supernova? type la, used to constrain the luminosity distance modulus, 



//(z): 



H(z) = 5 \og 10 (d L {z)) + 25; d L (z) = (l + z) 



Hz) 



-dz 



we evaluate statistical errors following the treatment in (1681 ) . Errors per each supernova 
are obtained adding in quadrature the uncertainty of the corrected apparent magnitudes 
due to the variation in the properties of SNe, <td, and the measurement uncertainty, 
a m . We neglect here the change in shape of the light curve of the Supernova? due to 
the variation in time of the gravitational constant (62), since by construction GR is 
the correct description of the Universe at late times. Systematic errors are estimated 
according to the prescription in (1691 ). as 

<j s (z) =A(1.7/,w)(l + z)/2.7, (12) 

where A is 0.05 for ground-based surveys, and 0.02 for space-based surveys, and z max is 
the depth of the survey. 

For the "2010" scenario we consider two cases: the Dark Energy Task Force Stage II 
for Supernova? (DETFII), as in Table 3.2 o f ([6 8), and the combination of the on-going 
ESSENCE (0), SNLS QH) and SDSSII Q surveys. We refer to such configuration 
as "ONGOING", and we assume a redshift distribution similar to the DETFII case 
(other than the local sample). For the "2020" scenario we consider a survey such as 
LSST or SNAP, corresponding to Stage IV of Ref. (1681 ) . 



BAOs. Baryon Acoustic Oscillations can be used to constrain separately the ex- 
pansion history, H {z)-iiom the line-of-sight clustering, and the comoving angular di- 
ameter distance, (1a(z) = adi,(z), if spectroscopic surveys are used. We follow (j?3 ) to 
forecast errors for both spectroscopic and photometric surveys; we always assume that 
systematic errors are below the statistical errors (]75l ). In the case of the BAO, the time 
line of "2010" is spread out over a few years. The first setup we consider is the SDSS 
LRG sample (1761 ); an improved version of the same survey (LRG BOSS) (1771 ). and the 
PAU-BAO survey (78). For the "2020" scenario we consider a survey like ADEPT (|79h . 



Weak lensing. Future weak lensing observations are used in order to constrain 
the power spectrum of the shear in multipole space, P K . In the Limber approximation 
it is written as (80): 



p„(Z) = ?js? g2{x) ■» ' 



where 



9(x) 



XII 



x X 



(13) 



(14) 
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and n(z) is the normalized source distribution, which we assume to be of the form 
n(z) = (z/zo) 2 x e^( 2 / 2 °^ 3/2 \ We compute the growth factor and matter transfer 
function using the analytical formulas of (81); we use the same approximation for the 
"true" JBD model matter power spectrum. We obtain the nonlinear matter power 
spectrum correction using the prescription of Peacock and Dodds (1821) . We however 
discard multipoles beyond I = 1500, in order to avoid errors coming from the uncertain 
nonlinear galaxy evolution and from baryonic physics (|83l ; [84j ). Errors on the shear 
power spectrum are obtained as 

(iLY 



6P K (l) 



p K (l) + 



+ Ode 



(15) 



/ (21 + l)/ sky v • • n / 
where 7; nt is the rms ellipticity per galaxy, assumed to be 0.16, h is the mean number of 
observed galaxies per square arcminute, and <7deg is the error coming from the photo-z 
errors in the determination of redshifts of sources. We estimate this error to be 1.6 x a z , 
where a z is the photo-z error of a give n survey, and we have used the analytic fitting 
formula for the shear variance (e.g. (I85I)). 

For the "2010" weak lensing scenario, we consider a configuration like DES (1861 ). and 
one like Pan-STARRS (1871 ). while for the "2020" scenario we assume a survey like DUNE 
Q0) or LSST (0). 



4-2. Results: forecasts for planned and future experiments 

We report specifics of the experimental setups considered and results for the statistical 
significance of detection of cujbd = 500 and 1000 in Tables EJ El HI For these two values 
of cujbd, we compute the quantities AdL, Ad a, AH and AP K , defined as in 

Aobs = obs(chain) — obs(true). (16) 

For each observable, "chain" refers to the value obtained using the cosmological 
parameters coming from the Monte Carlo chain in order to compute late time 
observables, and "true" is the observable for the fiducial JBD model. If gravity was 
GR, all the quantities above would be identically zero, because early and late time 
observables would be consistent with each other. Conversely, since we are assuming the 
the true Universe is a JBD, gravity modifications will manifest as a non-zero value for 
Ad L ,Ad A and AP K . 

We compute errors as follows, using Adi as an example. For all the models in the la 
vicinity of the chain, di (chain) will be larger than o^(true), because of the shift in the 
cosmological parameters as recovered by the chain. The error on di (chain) is evaluated 
as the difference between di for the best fit model and the minimum value of di for 
models within la from it; the error on o?L(true) is the experimental one around the true 
model, discussed in the previous Section. The total error on Ad L is obtained adding 
in quadrature errors on di (chain) and on ^(true). The statistical significance of a 
deviation of AdL from zero is given by the Ax 2 of AdL, evaluated using this total error. 

The value cujbd — 500 is found to be within reach of all the next-generation 
experiments (with a little delay for the BAO projects with respect to the other probes). 
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It is however interesting to compare how different setups can lead to different significance 
levels. For type la Supernova?, in the "2010" time range, we have considered two 
configurations, whose basic difference is the presence of ~ 350 more Supernovaae in 
the local sample for the DETFII case (|68l ) . Such drastic reduction of statistical error 
in the first bin would be expected to increase the statistical significance of detection. 
However, for our adopted estimate of Eq. (|12p for the systematic errors in ground-based 
surveys, having more than ~ 100 SNe per redshift bin of Az = 0.1 does not improve 
the signal-to-noise. This explain why the two surveys give similar performances. Both 
of them will be able to detect a JBD parameter ojjbd — 500 at the 3cr level. As for 
the "2020" scenario, a survey like LSST or SNAP will gain information both from the 
richness of the local sample and the deepness of the whole survey. In this case, we found 
that a JBD model with c^jbd = 500 can be distinguished from a GR, A CDM one with 
a significance of more than 7 a, and a value of ujbb = 1000 can be detected at the 3.3 
a level. 

For the BAO, we first consider constraints coming from both the comoving angular 
diameter distance, d^, which we anticipate to carry most of the signal-to-noise, and 
H(z). Errors are driven by the fraction of the total Universe volume covered by 
observations, so that a shallow (z max < 0.5) survey, such as the first setup of SDSS 
LRG, does not provide any significant detection of c<Jjbd = 500. Going up to z ~ 0.75, 
as in the SDSS LRG BOSS configuration, is enough to get a first detection at the 2.3a 
level, and the PAU-BAO project, which will reach a redshift of 0.9 and will cover 30% 
more square degrees, is competitive with the "2010" SNe probes described above. For 
comparison, in order to get a 3a detection of wjbd = 500, the shallow SDSS LRG survey 
should be able to observe over 30000 square degrees, more than 4 times its present sky 
coverage. In the "2020" experiments class, for a survey like ADEPT we forecast a 
detection of wjbd = 500 at the 7.3 a level, and a detection of cujbd = 1000 at 3.5 a. For 
both values, this result is at the same level, or slightly better, than the corresponding 
SNe experiments planned on analogous time scale. 

It is interesting to study what part of the signal-to-noise comes from information on H(z) 
along the line-of-sight, as allowed by spectroscopic surveys like the ones we consider. 
We find that the H(z) information is very useful for shallow surveys. In fact, with the 
usual notation, the difference 

AH = H (z) (chain) - H (z) (true) (17) 

is larger at low redshifts and almost vanishes at redshift larger than 1. As a result, the 
significance of this measurement is comparable with that of AcIa{z) for the LRG BOSS, 
or the PAU-BAO. On the other hand, for a survey like ADEPT, where information 
is collected at high redshift, measurements of H(z) do not substantially improve the 
performance of the experiment in detecting cujbd- 

For the weak lensing we forecast a significant detection of both cujbd = 500 and 
1000. In the case of cujbd = 500, we obtain a 4.1a significance of detection with an 
experiment like DES, and the number increases to 6.0 a for a survey like Pan-STARRS, 



Observational signatures of Jordan- Brans- Dicke theories of gravity 



14 



which should also be able to detect the value cjjbd = 1000 with a 3.6a significance. 
As for the "2020" scenario, a survey like DUNE or LSST will be able to reveal a value 
of ^jbd = 500 and cujbd = 1000 with a 6.8 and 3.8 a significance, respectively. The 
fact that Pan-STARRS and DUNE/LSST give similar results is due to the fact that 
for such large sample of galaxies, the main error on AP K (l) is due to the half-width of 
the la region within the best fit of the chain. The latter is only determined by the 
value of cjjbd, and does not depend on the survey. These results from lensing can be 
understood if we consider the known (approximate) dependence of the shear variance 
on the parameters f2 m , T, and a 8 (e.g. (1851) ). and study how the shear power spectrum 
varies with the shift of cosmological parameters. Combining the different contributions, 
we found that there is a strong (dlog P{k) / dlogVt m > 1) dependence on the shift in Q m , 
which is, as we have seen, the most sensitive parameter to the JBD field. Accordingly, 
the model reconstructed from the chain has ~ 10% more power at all scales than the 
"true" JBD model for cujbd = 500, and ~ 5% more power for c^jbd = 1000; these 
differences are larger than those found in all the other observables. 
However, as we already pointed out in the previous Section, we have only included the 
effect of some type of systematics. We have discarded modes above I = 1500, possibly 
eliminating sources of errors related to nonlinear evolution and baryonic physics, and 
we have accounted for the uncertainties in the determination of the source redshift; 
nonetheless, other types of systematic uncertainties, such as intrinsic alignment, PSF 
correction etc., may remain, possibly with amplitude comparable to the signal. For 
example, since the truly observable quantities are the correlation functions rather than 
the power spectrum P K , a further source of error is introduced in the mapping between 
the two (19 ll ). The fact that correlation functions cannot be observed over an infinite 
range of angles may bias the reconstruction of P K , which manifests as a fictitious 
oscillatory feature. The impact of such systematics is difficult to quantify extrapolating 
from the errorbars of the current, much smaller surveys; it may however contribute to 
degrade the efficiency of the weak lensing observations. 

For all the late-time observables we have considered, we show two relevant 
experimental configurations, one for the "2010" time line and wjbd = 500, and the 
other for "2020", u JBD = 1000, in Fig. H 



Table 2. Specifics of planned and future experiments for Supernova? type la, and 
statistical significance of forecast detection of wjbd = 500, 1000. The minimum redshift 
is assumed to be for all surveys, and we divide SNe in bins of width 0.1. 



Survey 


date av. 


max z 


SNe tot 


0"(^JBD = 500) 


a(cu JBD = 1000) 


DETFII 

ONGOING 

LSST/SNAP 


2010 
2010 
"2020" 


1.0 
0.9 
1.7 


1200 
850 
2500 


3.1 
2.9 
7.0 


1.4 
1.4 
3.3 








**** 



-4 + ' 




Figure 4. Comparison of for "true" JBD model (blue) and the extrapolation, assuming 
GR, of the la confidence region from the CMB chain (green). In left panels the 
JBD model has wjbd = 500, in the right panels wjbd = 1000. Top panel: distance 
modulus, rescaled to an empty Universe, for DETFII (left) and SNAP/LSST (right). 
Upper middle panel: log of the comoving angular distance, for the PAU-BAO (left) 
and ADEPT (right) surveys. Lower middle panel: Hubble parameter, for the PAU- 
BAO (left) and ADEPT (right). In the first case constraints from H(z) are in this 
case comparable to those from <1a{z), while the two sets of data points lie on top of 
each other for ADEPT'S redshifts, so that H (z) does not carry additional information. 
Bottom panel: Convergence power spectrum as a function of multipole, for DES (left) 
and DUNE/LSST (right) 
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Table 3. Specifics of planned and future experiments for BAO, and statistical 
significance of forecast detection of wjbd = 500, 1000. Where appropriate, we report 
in brackets the significance level without using the information from from H(z). 



Survey 


date av. 


min z 


max z 


deg 2 


o-(wjbd = 500) 


cj(cjjbd = 1000) 


SDSS LRG 

SDSS LRG (BOSS) 

PAU-BAO 

ADEPT 


2008 
2013 
2014 
"2020" 


0.14 
0.14 
0.1 
1.0 


0.47 
0.75 
0.9 
2.0 


7000 
7000 
10000 
30000 


2.3 (1.8) 
3.1 (2.6) 
7.3 (7.2) 


1.2 
3.5 (3.4) 







Table 4. Specifics of planned and future experiments for weak lensing, and statistical 
significance of forecast detection of wjbd = 500, 1000. 



Survey 


date av. 


med z 


deg 2 


n/arcmin 2 


photo-z err 


c(^jbd = 500) 


a{ojj BB = 1000) 


DES 


2009 


0.7 


5000 


10 


0.05(l+z) 


3.6 


2.1 


Pan-STARRS 


2012 


0.7 


30000 


5 


0.06(l+z) 


6.0 


3.6 


LSST/DUNE 


"2020" 


1.0 


20000 


100 


0.025(l+z) 


6.8 


3.8 



5. Conclusions 



We have studied the observational features of Jordan-Brans-Dicke (JBD) theories of 
gravity, as opposed to General Relativity (GR). Deviations from GR are described by a 
scalar field non-minimally coupled to the graviton; when the coupling parameter of the 
theory, cujbd, goes to infinity, the GR limit is recovered. We have proposed a method 
which would allow one to reveal failures of GR without assuming previous knowledge of 
the true theory of gravity. 

We have first exploited the well-known method of assuming the same expansion 
history for a JBD and a GR model, and analyzing the associated perturbation growth, 
which will depend on the underlying theory of gravity (1231 ; |24J; |25|; |26j; |27j; |28|; |29|; l30l; l3ll). 
We showed that this method is not really applicable in scenarios like this, where 
gravity modifications take place at early times. This suggested the idea of looking 
for inconsistencies between "early" and "late" time observables as opposed to expansion 
history and perturbation growth. 

The method we propose is general and can be used to discriminate between any 
two theories who give similar predictions at some epoch, and diverge at some other 
epoch. It only relies on the fact that a correct extrapolation of observations made at 
early times to late times, or vice- versa, requires the knowledge of the theory of gravity. 
Even if gravity is GR, it is still possible that the extrapolation gives wrong results if 
the matter-energy content of the Universe is mistaken, such as in Early ( 1381 : l39l : |40j; l4ll ) 



versus Ordinary (1351 ; l36l ; 1371 ) Quintessence models. 
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We have used the primary CMB as the early-time observable and Supernovas, Baryon 
Acoustic Oscillations and Weak Lensing as the late time observables. Further gain in 
significance of the results could be obtained considering more observables. The "early- 
times" class is, of course, the more troublesome as for direct observations; in this respect, 
future surveys of the 21cm hydrogen transition line from the high-redshift IGM (j92j) 
could provide relevant improvements. 

We translated our results into limits on the coupling parameter of the theory, cujbd- 
We showed that next-generation experiments will be able to improve substantially the 
current limit cujbd > 120 (1221 ); detection of values as large as 500 and 1000, respectively, 
are within reach of the "2010" and "2020" experiments, for all the observables we 
considered (SNe, BAOs, Weak Lensing). 

Limits on the JBD parameter can be thought of in terms of an effective 
parametrization of deviation from GR and in particular can be interpreted as limits 
on the variation of the gravitational constant, G. In fact, we have seen that the 
evolution of the effective gravitational constant is driven by the evolution of the JBD 
field, so that these observations constrain the time variation of fundamental constants 
(1201 : 1931 ) . Constraints of the type cjjbd > 500, 1000 would correspond to a variation 
(|G rec - G \)/G < 1.38% and 0.69% respectively. 

We also note that all the observables used in this paper are among those indicated 
as primary science goals by the recent reports (]94j ; l68l ) . 

To conclude, let us recall that limits on cujbd on cosmological scales, as in (1221 ) and 
previously in (1531 ). were obtained starting from the assumption of JBD as the correct 
theory of gravity. We also aim to compare constraints obtained with such approach and 
with the one we used in this paper, on the basis of the same presently available data 
sets (JotI ). Such analysis will help to clarify which is the best approach to pursue in the 
quest for the true theory of gravity. 
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6. Appendix A: Behavior of the matter density contrast in GR and JBD 

We have seen that our formulation of JBD and GR differ mostly at early time and that 
therefore late-time observables, such as the density contrast 5 m , are not sensitive to 
ujbd- We justify such statement here in more detail. 

The redshift dependence of S m (z) in GR can be obtained from its evolution equation: 

8 m {z) + 2H(z)5 m (z) = AixG P (z)5 m (z). (18) 

H(z) determines p(z) exactly through the Friedmann equation, so that the coefficients 
of the above equation are unambiguously defined. 

On the other hand, we have learned that the modifications to GR that we are considering 
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induce a time dependence of the gravitational constant, so that the above equation 
becomes (95) [jj 

8 m {z) + 2H{z)6 m {z) = 4nG cS (z)p(z)6 m (z). (19) 

The solution for 5 m (z) in this case, even assuming the same H(z), will be different from 
GR for two reasons: first, the mapping between H(z) and p(z) will change, as seen in 
the previous section, and second, the source term of such equation will also change, in 
response to the time variation of G e g(z). 

However, how well can we expect this method to do for JBD models? The function 
G e ff(z) in this case is l/(87rF): this implies, given the evolution of 1/F from Fig. [fl 
that the time dependence of H(z),p(z) and G e s(z) will differ significantly from the 
ordinary gravity case only at high redshift. At redshift, say, z < 1, the coefficients of 
the differential equations ( fl8|) and ( fl9i) will be similar, and so will be the solution S m . 
To illustrate all this, we have developed a method to reproduce exactly the expansion 
history of a given JBD model, which we will call H(z) ie f, in an Ordinary Gravity case. 
This cannot be done by simply changing the cosmological parameters: in a flat A CDM 
model H(z) only depends on /i and Q m (and fl r for z » 100), 

H GR (z) = H 1 , (20) 

ln r (l + z) 4 + Q m (l + z) 3 + Q A 



But in the JBD case there is no such simple functional form of the Hubble factor; the 
complicated evolution of H Te f(z), given by Eq. [TU], cannot be described as a sum of the 
three power-laws in redshift which appear in the equation above. 

We consider a GR, three-fluid model with ordinary matter and radiation components 
plus a perfect fluid with equation of state u>jbd(z). We will refer to such component as 
the "Jordan-Brans-Dicke fluid", and other than the value of u>jbd we will impose the 
same cosmological parameters as the "true" JBD model. This way, we expect to isolate 
the effect of the JBD field at early times and to recover the standard Cosmological 
Constant term at late times. For this model, the expression for H (z) is identical to Eq. 
(120]) but with the substitution — ► ^jbd(^) where 



^jbd(^) = ^jbd(^ = 0) exp 



f° 

-3 J (1 + w JBD (z))dz 



(21) 



^jbd(-2 = 0) is the present value of the JBD field density, and can be fixed as 1 — Q m 
requiring geometrical flatness. 

With a generic u>jbd( z ), any function H(z) can be exactly reproduced in this form. The 
relative weight of the "JBD fluid" density measures the contribution of the JBD field 
to the total energy density, and therefore to the expansion history. It can be evaluated 
as (pjbd — Ptot)/Ptot; and is shown in the left panel of Fig. [5] for the expansion history 
of a model with cjjbd = 100. On the right panel we show the equation of state of the 
fluid, u>jbd; which, is helpful for the interpretation of physical effects of the JBD field. 

§ We are indeed neglecting anisotropic stress, which would introduce a non-trivial dependence of the 
solution on the wavenumber k. However, we have numerically checked that the difference of the two 
gravitational potentials is very small in such models, so that the error in this approximation is < 10 -4 . 
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Figure 5. Left panel: Relative weight of the hypothetical fluid component in a 
minimally coupled model with respect to a JBD model with cjjbd = 100 and coinciding 
expansion history. Right panel: Equation of state of such fluid component. 



At very early times the relative contribution of such component is five to ten percent 
for this value of cujbd, and its equation of state is positive and in the range between 0.2 
and 0.3. Such fluid may be interpreted as a mixture of matter and radiation. 
At the CMB redshift, z ~ 1080, the equation of state of the JBD fluid is close to 
zero, mimicking an additional matter component, which is indeed expected from the 
behavior of the multiplicative term 1/F in the density equation. We have seen that at 
early times the gravitational coupling in the JBD model is stronger than in GR; in order 
to reproduce this feature in the Ordinary Gravity case, one would need to enhance the 
total matter density. 

At late times the fluid behaves, as anticipated, as a Cosmological Constant component, 
and its equation of state at z < 1 is practically indistinguishable from —1 . Therefore, 
at late times we cannot expect such fluid to contribute to the evolution of S m . To 
illustrate this, in Fig. [6] we plot the power spectrum of matter, evaluated at the 
present time, for the two models, and the ratio of their perturbation growth factors. 
Although for wavenumbers k > 0.1 Mpc -1 there is a relevant difference between the two 
power spectra, it has to be attributed to the difference in the transfer function. The 
field behaves as an additional matter component, therefore the horizon size at matter- 
radiation equality changes, and so does the transfer function on smaller scales. The 
redshift evolution of the growth factor g(z) in the two cases only shows differences of 
< 4%, and, as expected, mainly at redshift higher than those probed by large-scale 
structure or lensing surveys. Furthermore, let us stress that the value of cujbd we are 
conside ring for illustrative purposes is fairly large and already ruled out by cosmological 
probes (1221 ) . It is indeed true that while techniques used so far to trace S m , such as galaxy 
surveys, weak lensing, cluster counts have only been used at z < 1, different methods may 
be used to probe much larger redshifts. The proposed technique of tracing the neutral 
Hydrogen using its 21 cm transition is challenging (e.g. 092!)), but very promising and 
may provide accurate maps of the high redshift large scale dark matter distribution. 
However, for the time being, the detection of deviations from GR for a JBD model 
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Figure 6. Left panel: Linear matter power spectra of the JBD model with wjbd = 
100 (solid line) and of a GR, ACDM model with the same expansion history and 
cosmological parameters. Right panel: Ratio of the perturbation growth factor g(z) 
for the same two models. 



through inconsistencies in the expansion history and perturbations growth are limited 
by the poor measurements of 5 m at early times. 

7. Appendix B: Code comparison 

The DEfast code, used to produce the reference JBD spectra, and the CAMB code used 



by the COSMOMC program are known to be in agreement within the 1% level (1961 ). 
However, the shift in the cosmological parameters we are looking for are of comparable 
size. In order to properly calibrate outputs, we run several chains, for different values of 
the cosmological parameters, for GR ACDM models, using input spectra from DEfast. 
Only if all the values of the cosmological parameters are recovered correctly one can 
be sure that the observed shift is due to the JBD nature of gravity rather than to a 
numerical effect. We report in the second column of Table [5] the result of one such 
chains, for the same reference model quoted in the main text. It can be seen that, for 
all the parameters, the values resulting from the fit lie within the 1 — a vicinity of the 
true ones. However, we noticed a slight positive shift in the reconstructed value of n s , 
of the order of 0.4%, which appeared to be of systematic nature. We thus compared 
the output of the DEfast code and the CAMB code for a sample of 520 models of the 
chain mentioned above, indeed finding some residual numerical difference. In particular, 
fixing the primordial amplitude normalization scale at k ~ 0.05 Mpc -1 , corresponding to 
l ~ 700, the DEfast code predicts a slightly lower power in temperature, with respect to 
CAMB, at lower multipoles, and slightly more power at higher multipoles. We however 
discovered that such difference is largely independent of the cosmological model: we 
averaged the ratio of the power spectra given by CAMB and DEfast over N = 10, 100 
and 520 models and found extremely similar curves, as plotted in Fig. [7J Let us notice 
that the shape of such difference is very close to what would be caused by a shift in the 
value of the spectral index of the order of the one we observe; we plot the corresponding 



REFERENCES 



25 



1.010 




Figure 7. Left panel: Ratio between the temperature power spectra as a function of 
multipole, for 10 (dotted line), 100 (dashed line), and 520 (solid line) sample models. 
Right panel: distribution of the x 2 due to the disagreement between DEfast and 
CAMB, for the temperature (dashed line) and EE polarization (solid line) spectra, 
as a function of multipole. 



to be compared with our empirical correction curve, in the same Figure (smooth solid 
line) . 

We concluded that the ratio of the power spectra from the two codes, averaged over a 
large number of models, can be used as an effective correction to the DEfast temperature 
power spectrum. Results of a chain run with the same cosmological parameters and the 
corrected Cf T are shown in the third column of Table [5j we see that the shift in the fitted 
value of n s is now negligible. This correction was used throughout all our numerical 
analysis. 

However, even after this correction, some small differences in the output from DEfast and 
CAMB remained. Even if not harmful for our procedure of reconstructing cosmological 
parameters, they are responsible for the fact that the value \ 2 f° r the reconstructed 
model does not approach zero as we recover GR. In fact, the fit to a GR ACDM model 
whose input power spectra are generated with DEfast has a value of \ 2 °f 9.06; this is 
due to the fact that COSMOMC uses CAMB in order to generate the CMB spectra. 
The same fitting procedure, if done using input spectra generated by CAMB, would 
yield a much smaller y 2 of 0.26. We checked the distribution in the multipole space of 
the x 2 ? finding that it is in large part due to slight differences in the EE polarization 
spectrum at low 1. The total value of the x 2 for this spectrum alone would be 5.15, while 
the temperature power spectrum would only have 0.91 (and their cross-correlation, of 
course, is responsible for the missing 3.02). We show the distribution of the \ 2 i f° r the 
TT and EE spectra in the right panel of Fig. [7J 



effect on the power spectrum: 




(22) 
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Table 5. Best fit and la confidence levels of cosmological parameters, before and after 
correction. 





"true" 


uncorrected 


corrected 


<^CDM 

n s 

T 

A 


0.022 
0.1232 
0.95 
0.09 
3.1355 
0.72 
0.28 


n D99 0.0222 
u - uzz 0.0218 
n lOQc; 0.1263 
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n 0990.0222 
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